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Abstract. In classifying the ensemble of powerful extragalac- 
tic radio sources, considerable evidence has accumulated that 
radio galaxies and quasars are orientation-dependent mani- 
festations of the same parent population: massive spheroidal 
galaxies containing correspondingly massive black holes. One 
of the key factors in establishing this unification has been the 
signature of a hidden quasar detected in some radio galaxies in 
polarized light. The obscuration of our direct view of the ac- 
tive nucleus usually, but not necessarily exclusively, by a thick 
nuclear disk or torus can act conveniently as a 'natural corono- 
graph' that allows a much clearer view of the host of a radio 
galaxy than of a quasar. 

In this study, we exploit the opportunity to eliminate the 
quasar glare by performing sensitive spectropolarimetry with 
the Keck II telescope of a sample of radio galaxies with red- 
shifts around 2.5. This represents the epoch when quasars were 
many times more common that they are now and is likely to 
be the period during which their host galaxies were being as- 
sembled into what become the most massive galaxies in the 
Universe today. We show that dust-reflected quasar light gen- 
erally dominates the restframe ultraviolet continuum of these 
sources and that a highly clumped scattering medium results 
in almost grey scattering of the active galactic nucleus pho- 
tons. The observations, however, do not exclude a substantial 
star formation rate averaged over a Gyr of evolution. The sub- 
mm reradiation from the scattering dust is likely to represent 
only a small fraction (~ 10%) of the total far infrared lumi- 
nosity. An analysis of the emission lines excited in the inter- 
stellar medium of the host galaxy by the hard quasar radiation 
field reveals evidence of a dramatic chemical evolution within 
the spheroid during this epoch. Secondary nitrogen production 
in intermediate mass stars produces a characteristic signature 
in the Nv/Civ and Nv/Hell line ratios which has been seen 
previously in the broad line region of quasars at similar red- 
shifts. We find intriguing correlations between the strengths of 
the Lya and Nv emission lines and the degree of ultraviolet 
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continuum polarization which may represent the dispersal of 
dust associated with the chemical enrichment of the spheroid. 



Key words: Galaxies: evolution; Galaxies: abundances; Galax- 
ies: active; Galaxies: starburst; Polarization; Scattering 

1. Introduction 

Radio galaxies are thought to be the hosts of radio quasars, 
oriented such that the direct view to their nuclei is obscured 
by an optically very thick nuclear structure (Antonucci 1993 
and references therein). This fortuitous occultation of the ac- 
tive galactic nucleus (AGN), which would otherwise outshine 
an galaxy by a factor of up to a thousand or so, acts as a 
convenient 'natural coronograph', allowing the detailed study 
of many properties of the surrounding stars and interstellar 
medium (ISM). Our current belief is that the hosts of powerful 
radio sources in the distant Universe are destined to become the 
giant ellipticals of today: the most massive galactic systems we 
know (McLure et al. 1999). While some may have commenced 
their formation at very high redshift, it is clear that the pro- 
cess of assembly is very active at z ~ 2.5 (Pentericci et al. 
1999, 2000 ; van Breugel et al. 1998). This corresponds to the 
epoch when luminous quasars appear to have had their maxi- 
mum space density (see eg. Pei 1995; Shaver et al. 1996; Fan 
et al. 2000). Indeed, a major goal of our programme is to un- 
derstand any causal relationship between the formation of the 
massive black hole and that of the galaxy within which it re- 
sides. 

Although we may not see it directly, the presence of the 
huge luminosity radiated by the quasar makes itself felt in 
many ways in the surrounding galaxy and understanding those 
processes is essential if we are to deduce the basic properties 
and evolutionary state of the host. The radiation field ionizes a 
substantial fraction of the gaseous interstellar medium (ISM), 
producing emission line halos extending sometimes hundreds 
of kpc (van Ojik et al. 1997a) and both destroys and acceler- 
ates dust (de Young 1998). The ultraviolet (UV) and optical 
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components of the quasar continuum and broad emission lines 
can now, thanks to large telescope polarimetry, be seen scat- 
tered from dust and/or electrons far from the nucleus (see eg. 
di Serego Alighieri et al. 1996a; Dey et al. 1996; Knopp and 
Chambers 1997b; Cimatti et al. 1998b). Both of these repro- 
cessing mechanisms, fluorescence and scattering, contribute to 
the alignment between the observed optical and radio struc- 
tures seen in radio galaxies beyond a redshift of 0.7 or so 
(Chambers et al. 1987; McCarthy et al. 1987). Although at op- 
tical wavelengths the effects of scattering are seen most clearly 
at higher redshifts, where the dilution by starlight below the 
4000 A break is small, the process can be seen to operate in 
low redshift radio galaxies when they are observed with large 
telescopes (Cohen et al. 1999). 

Above 4000 A in the restframe, the evolved stellar popu- 
lation contributes substantially to the observed continuum and 
produces an infrared (K-band) Hubble diagram with a remark- 
ably small dispersion (McCarthy 1993, but see references for 
the higher redshift sources, eg. Eales et al. 1997; van Breugel 
et al. 1998). Below 4000 A, we might expect to see directly 
a young stellar population in a starburst arising as a precur- 
sor to the AGN or, alternatively, as a result of its influence 
on the surrounding medium. Objects with a redshift around 
z ~ 2.5 give optical spectrographs access to the restframe spec- 
tral band which includes the strong UV emission lines from 
Lya through Cm] and beyond and also to that part of the con- 
tinuum which displays detectable spectral features from OB 
stars as well as absorption lines from the ISM. In addition, the 
spectrum covers the region where Galactic-type dust displays 
the 2200 A signature and the steep extinction increase below 
1600 A. 

Our programme of Keck observations of z ~ 2.5 radio 
galaxies is designed to obtain high quality spectrophotome- 
try of the spectral features mentioned above together with po- 
larimetry having sufficient signal-to-noise to detect a contin- 
uum polarization as low as one or two percent. The polarimetry 
is used to separate the directly received light from the scattered 
radiation while the low-noise line and continuum spectroscopy 
can be used to study many properties of the host galaxy, es- 
pecially the physical and dynamical state and the chemical 
composition of the ISM. Data from the first two sources ob- 
served, 4C+23.56 and 4C— 00.54, are presented in Cimatti et 



al. (1998a) (hereafter paper I). Here, we report the extension 
of the sample to nine sources. In this paper, the analysis is re- 
stricted to the integrated properties defined by the continuum 
extent. An analysis of the spatial variations will come in a later 
publication (Villar-Martm et al. in prep.). Throughout this pa- 
per we assume Ho = 50 km s _1 Mpc -1 and go = 0.1. 

2. Observations and data reduction 

The list of objects and the journal of observations are given in 
table [Ij All sources were selected from the ultra-steep spectrum 
(USS) radio galaxy survey (see eg. Rottgering et al. 1995b) 
with redshift greater than 2.2 in order to observe the Lya line. 
While the statistical properties of this subset cannot be defined 
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Table 1. Journal of observations. (1) IAU object name; (2) 4C 
name (3) Redshift; (4) Total exposure time; (5) slit position 
angle; (6) Run a: July 05-06 1997 (data published in Cimatti et 
al. (1998)) , b: Dec. 24-25-26 1997, c: May 26-27 1998 



with precision, it is effectively unbiased with respect to both ra- 
dio flux and optical magnitude for the identified sources (Car- 
illi et al. 1997). The R magnitudes range from 21 to 23 (see 
table ^, column 8). The source 4C+03.24, which has a redshift 
of 3.6, was originally included as an attempt to investigate the 
dependence of properties on epoch but, due to the faintness of 
many of the more distant objects, it was subsequently decided 
to restrict the sample to z < 3. 

Observations were made using the Low Resolution Imag- 
ing Spectrometer (LRIS, Oke et al. 1995) with its polarime- 
ter (Goodrich et al. 1995b) at the Keck II 10m telescope from 
July 1997 to May 1998 under subarcsecond seeing conditions 
(seeing ranging from 0.5" to 1"). The LRIS detector is a Tek 
2048 2 CCD with 24/im pixels which correspond to a scale of 
0.214"pixel _1 . We used a 300 line mm -1 grating and 1" wide 
slit which provide a dispersion of 2.4 A pixel -1 and an ef- 
fective resolution of ^10 A FWHM. The spectral range is 
A 6s ~ 3900 — 9000 A. Observations of each object were di- 
vided into several sets (between 2 and 5) of 4 exposures of ap- 
proximately 30 minutes with half-wave plate position angle set 
successively to 0°, 45°, 22.5° and 67.5°. The slit was always 
oriented along the radio axis (Rottgering et al. 1994, Carilli et 
al. 1997a). 

After debiasing and flatfielding, 2D-spectra were rectified 
for geometrical distortion and wavelength calibrated using Hg- 
Kr arc lamp spectra. Spectra were cleaned for cosmic rays 
hits and then extracted in apertures given in table || in order 
to include all the extended continuum flux. After performing 
a precise wavelength registration of all eight spectra in each 
set (ordinary and extraordinary rays for each half-wave plate 
orientation) using night-sky emission lines to determine final 
wavelength zero-points, data were binned and then combined 
as described in Cohen et al. (1995, appendix A) to form I, Q 
and U Stokes parameters. Unbiased values for the fractional 
polarization were estimated with the best estimator given by 
Simmons and Stewart (1985). Statistical confidence intervals 



on fractional polarization and polarization angle were deter- 
mined using a Monte-Carlo method taking into account the 
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background polarization and the detector noise (details of the 
method in Vernet 2000). Polarized (VICygl2, BD+332642) 
and unpolarized (HZ44, BD+284211) standard stars were ob- 
served in order to check and calibrate the polarimeter. Flux cal- 
ibration was done with the spectrophotometric standard stars: 
GD248, G191-B2B, Feige 22, Feige 66, Feige 67. Since sev- 
eral nights were not photometric, the spectra have been scaled 
to published HST magnitudes when available (Pentericci et al. 
1999 , see column 8 in table ^) without any aperture correction. 
Spectra were corrected for Galactic extinction using Burstein et 
al.( [l982]) maps and adopting the extinction curve of Cardelli et 
al.p89l 



A composite spectrum of all nine radio galaxies was also 
computed. The result is shown in fig. [|. We first normalized 
the spectra in a line free region between 2000 and 2100 A and 
then combined them using a simple average. 



3. Results 

3.1. The continuum 

The spectropolarimetry results are presented in fig. [j], On each 
panel is displayed, from top to bottom: the total flux spectrum 
at two different scales, the fractional polarization and the posi- 
tion angle of the electric vector. The main results are summa- 
rized in table ||. The aperture used for the extraction matches 
the continuum spatial extent (table ||, column 2). The contin- 
uum is clearly resolved in all objects extending from 4" up to 
8" which corresponds to a linear scale at z ~ 2.5 of 45 to 90 
kpc. It shows a variety of morphologies such as double narrow 
sources (4C+23.56), very extended emission (4C +48.48) and 
single peaked sources (0943—242). This paper does not dis- 
cuss further the spatial structures but concentrates solely on the 
properties integrated over the regions where the continuum is 
detected. 

The first important result is that all high reshift radio galax- 
ies (HzRGs) in this sample except 0731+438 show high con- 
tinuum fractional polarization ranging from ~6% up to ~20% 
measured in a large bin just longward of Lya (see table ||, col- 
umn 3). This window, between Nv and the Olv]/SilV blend 
(~ 1250 — 1400 A in the restframe), contains several weak, 
low ionization absorption lines but is the broadest gap between 
strong emission lines in a region of the spectrum where the 
instrumental sensitivity is high for all sources. For 0211 — 122 
and 0943—242, we clearly see a rise of the continuum polar- 
ization toward the blue. In 0211 — 122, the fractional polar- 
ization rises from P = 15.3 ± 1.1% shortward of Cm] to 
P = 19.3 ± 1.2% longward of Lya and for 0943-242 it rises 
from P = 4.1 ± 0.6% to P = 6.6 ± 0.9% berween the same 
rest wavelengths (note here that error bars are symetric because 
data were binned in order to reach signal-to-noise P/a > 6). 
A similar trend is also seen in 4C+23.56 and marginally in 
4C— 00.54 (see paper I). The signal-to-noise of the polarization 
measurements is generally not high enough to detect any sig- 
nificant variation of fractional polarization with wavelength in 
other galaxies in our sample. For 0731+438, we can only give 
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Fig. 2. Histograms of the difference between PA ra d and PA po / 
(PA of the E- vector, light) and between PA^v and PA po / (dark). 



an upper limit of about 2.4% at the la level to the fractional po- 
larization just longward of the Lya line. However, a significant 
polarization is measured at longer wavelength (P = 3.1+0.7% 
shortward of Cm]). 

The position angle (PA) of the electric vector (PA po /, table 
||, column 4) is generally constant with wavelength within the 
measurement errors. We show in fig. || the difference between 
the radio axis position angle (PA ra d) and the polarization an- 
gle (PA of the E-vector, light) and the difference between the 
UV extension position angle (PA^y*) and the polarization an- 
gle (dark). Averaging over the nine sources in our sample we 
find PA ra d— PApo/ — 69°. For the 6 sources for which we have 
PAjjv measurements the average PAyy— PA po / is 84°. This is 
consistent with the PA of the electric vector being generally 
perpendicular to the UV continuum extension but less closely 
correlated to the large-scale radio structure (this effect was al- 
ready noted by di Serego Alighieri et al. 1993, Cimatti et al. 
1993,1994 and Hines & Wills 1993). 

One very striking result is that all objects in this sample 
show a remarkably similar continuum shape with a dip around 
2200 A and a shallow peak in F\ between Lya and 1500 A. 
This similarity of shape and in particular the 2200 A dip can- 
not be an artifact of the response calibration since objects in 
the sample cover a relatively large range in redshift. For in- 
stance, 0943-242 at z = 2.922 (1000 < X rest < 2250) does 
not show any dip between observed wavelength 7000 A and 
8000 A. Although spectra are quite noisy at long wavelength, 
the continuum rise at A resi > 2200A is well detected using 
just regions between strong sky lines (see in between shaded 
areas showing region of strong sky emission in fig. |l|). The de- 
crease below 1500 A is also a real continuum feature since 
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Fig. 1. Spectral and polarization properties. In each panel, from top to bottom: the observed total flux spectrum in units of 
10~ 17 erg s _1 cm~ 2 A -1 plotted at two different scales, the first to show strong emission lines and the second to show the 
continuum, the fractional polarization in % and the position angle of the electric vector (measured N through E). Filled circles 
and crosses respectively indicate continuum and narrow emission lines (including their underlying continuum). Shaded areas 
indicate regions of strong sky emission. Left : 4C +40.36; Right : 021 1 — 122 
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Fig. 1. Continued. 4C+03.24 
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Table 2. Continuum properties. (1): Object name; (2): Aperture used for spectral extraction (the slit width was 1" except for 
4C+23.56 and 4C— 00.54 for which the slit was 1.5" wide) ; (3): Fractional polarization computed in a large bin between Nv 
and SilV lines (~ 1250— 1400 A) with 1 a confidence interval (bins are not necessarily identical to the ones shown on fig. [I]) ; (4): 
Electric vector position angle with 1 a confidence interval ; (5) Position angle of the UV continuum extension as measured on 
HST images (Pentericci et al. 1999); (6) Slope /3 between 1500 A and 2000 A (with F\ oc A' 3 ) ; (7) Logarithm of monochromatic 
luminosity at A rest «1500 A in erg s^A^ 1 ; (8) R band magnitude from Carilli et al. (1997) and Rottgering et al. (1995); (9) 
HST magnitude from Pentericci et al. (1999). 



it starts before the onset of the Lyman forest absorption below 
Lya. Continuum slopes between 1500 A and 2000 A are listed 
in table |2|, column 6. 



3.2. Emission lines 

Fluxes of all detected emission lines are given in table ^. The 
quality of the data allow us to detect a number of weak emis- 
sion lines such as SiiiA1309, Niv]A1488 and Nm]A1750. The 
spatial extent of strong emission line halos (Lya, Nv, Civ and 
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Hell) ranges from about 45 kpc in 0943—242 up to more than 
100 kpc in 0211-122 and 0828+193 (see also van Ojik et al. 
1996). These strong emission lines are usually more extended 
than the aperture used for the extraction and show very com- 
plex spatial and velocity structure with high velocity ampli- 
tudes up to ^2000 km s -1 . Narrow emission lines have typical 
width of ^800-1500 km s -1 as measured on the HellA1640 
line (see table which is less affected by possible underlying 
broad emission (since it is usually very weak in quasars) or ab- 
sorption. Broad lines are clearly detected in 0211 — 122 (broad 
Civ FWHM - 6500 km s ~\ W\ ~ 57 A and broad Cm] 
FWHM ~ 9000 km s _1 , W\ ~ 50 A) and 4C+23.56 (broad 
Civ FWHM ~ 4800 km s~\ W x ~ 68 A and broad Cm] 
FWHM ~ 6300 km s" 1 , W x ~ 76 A). Broad wings are also 
marginally detected in 4C +40.36 and 4C -00.54. 

Since narrow emission line profiles can be quite poorly 
represented by gaussians due to the presence of several kine- 
matic components, we did not attempt to perform gaussian 
fitting. Fluxes were measured from simple line integrals. The 
main source of error in these measurements comes from uncer- 
tain continuum determination. We estimated the uncertainty by 
making limiting high and low estimates of the continuum level. 
This method does not provide rigorously defined measurement 
errors but does give a conservative estimate of the range in flux 
allowed by the data for each emission line. 

One particular issue in measuring narrow emission line 
fluxes is the contamination by broad components, especially 
for the Lya, Nv, Civ and Cm] lines. In the two cases where 
broad lines are most prominent, (4C+23 .56 and 02 1 1 — 1 22) and 
most likely to significantly affect our results, we made mea- 
surements of the narrow emission lines using several different 
strategies including multiple-gaussian fitting, integration from 
an elevated continuum and measurement after subtraction of a 
scaled quasar template. The small dispersion in the resulting 
values is included in our conservative estimate of the uncer- 
tainties. 

We also measured the polarization of the strongest narrow 
emission lines including the underlying continuum (see bins 
marked with crosses in fig. [j]). They show, in general, a signif- 
icantly lower polarization than neighbouring pure continuum 
bins. In all cases the continuum subtracted polarization in the 
narrow Lya line is less than 2%. 

3.3. Absorption lines 

Spectra generally show several absorption features that can be 
identified as low ionization zero-volt or excited fine-structure 
interstellar lines such as Sill AA 1260,1265, Ol AA 1302,1305 
and Sin AA 1304,1309, Cn AA 1334,1336, Sin AA 1527,1533. 
There is a weak tendency for high polarization objects to show 
these lines in emission and low polarization objects to show 
them in absorption. Equivalent width measurements and iden- 
tifications of these features are given in table [7| Fig. ^illustrates 
the complex behavior of these lines with the two extreme cases 
observed in our sample. While low ionization lines are only 
seen in absorption in 4C+23.56 (bottom spectrum in fig. ||), 
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Table 3. Equivalent width and instrumental profile (10 A) 
corrected observed FWHM measurements of the Hell A 1640 
emission line. FWHM measurements were done using a sin- 
gle gaussian fit in all cases without any attempt to deblend the 
different kinematic components. 
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Fig. 3. Low ionization absorption and emission features for two 
extreme cases. While low ionization lines are only seen in ab- 
sorption in 4C+23.56 (bottom spectrum) these lines only show 
up in emission in 0943—242 (top spectrum). The spectrum in 
the middle is our HzRG composite spectrum. 



these lines only show up in emission in 0943—242 (top spec- 
trum in fig. ||). 

In the case of 0731+438 and 4C-00.54 (see paper I), an 
absorption is detected close to the expected position of the 
Cm AA 1426,1428 photospheric line but we consider this iden- 
tification as very uncertain since no other absorption lines com- 
patible with pure photospheric lines are detected. See sect. 
4.1.3 for a more detailed discussion on photospheric absorp- 
tion lines. 



3.4. Notes on individual objects 

- 0828+193: this radio galaxy was thought to have two sub- 
components separated by 3" (Oil to the SE and OI2 to the 
NW, see Knopp and Chambers 1997). Our deep spectrum 
reveals that Oil is in fact a M5V type star. The 3a up- 
per limit to the polarization of this star between 5600 and 
8300 A is 0.85 %. 
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Fig. 4. Spectrum of intervening galaxies (upper spectrum), left : Hll galaxy 2 " southwest of 4C+40.36 with template galaxy 
NGC 6052 overlayed; right : Hll galaxy 2.5 " southwest of 4C+48.48 with template galaxy NGC 6764 reddened with Eb-v = 0.1 
The main emission lines are labeled. In both cases, spectra have been cleaned of residual extended emission lines from the radio 
galaxy and smoothed with a 5 pixel boxcar. Template galaxies were obtained from Schmitt et al. (1997) and have been shifted 
downward by 10 flux units for clarity. Shaded areas show regions of strong sky emission. 



- 4C +40.36: Chambers et al. (1996b) note that there is a 
small extension 2" to the southeast of the eastern lobe 
which is detected in all bands and has very similar colours 
to 4C+40.36. They identify it as a good candidate for a 
companion or a subcomponent (see their fig. 7). This is in 
fact a foreground Hll galaxy at z = 0.404. The spectrum 
is displayed in fig. ^a with an overlay of the Hll galaxy 
NGC 6052 from Schmitt et al. (1997). 

- 4C+48.48: we find an emission line object 2.5 " to the 
southeast of 4C+48.48 that has a significantly redder con- 
tinuum than the radio galaxy. The measured polarization 
between 5600 and 8000 A is compatible with 0% with a 
3(7 upper limit of 4.5%. A strong emission line is clearly de- 
tected at A hs = 6277.5 A. We also identify another weaker 
emission line at A /, s = 8184.7 A. These lines do not corre- 
spond to any expected strong emission at the redshift of the 
radio galaxy. Identifying these lines as [Oll]A3727 and Hp, 
this intervening object can be classified as an Hll galaxy at 
z = 0.684. The spectrum is displayed in fig. |]b with an 
overlay of the Hll galaxy NGC 6764 from Schmitt et al. 
(1997). 

Note that, due to the presence of intervening objects 
along the line of sight, observed properties of 4C+40.36 and 
4C +48.48 might be significantly affected by lensing. 

3.5. Additional data 

We also collected information on three high redshift sources 
that have comparable quality spectropolarimetry data from the 
literature: two radio galaxies (4C+41.17 and 3C 256) and the 
ultraluminous infrared galaxy FSC 10214+4724. The main 
properties of these objects that are used in the discussion are 
given in table 



Six sources in our sample have been observed with 
HST/WFPC2. These data are published in Pentericci et al. 
(1999). They were mainly used in this work to compare the 
polarization angle to the direction of the extended ultraviolet 
emission (see table and to obtain an absolute flux calibration 
of our spectra. 

Radio properties of all objects in our sample can be found 



in Rtittgering et al. (1994) and Carilli et al. (1997) 



4. Discussion 

The discussion is based on the integrated properties of the 
sources. The spatial properties and kinematic details are ad- 
dressed in a separate paper (Villar-Martm et al. in prep.). The 
aim is to identify the origins of the various radiated continuum 
and line components using polarimetry as a powerful discrimi- 
nent between direct and scattered radiation. 



4.1. The nature of the continuum 

The question we address in this section is the origin of HzRG 
ultraviolet continuum: what is the dominant source of UV con- 
tinuum and what is the relative contribution of the different 
possible components? The high fractional polarization together 
with the presence of broad lines clearly suggest that scatter- 
ing of the radiation from the central active nucleus plays an 
important role, and probably a dominant one in the most po- 
larized objects. The universality of the shape of the spectral 
energy distribution is difficult to reconcile with the dominance 
of a young stellar population in this wavelength range since 
starbursts show strong colour evolution and a range of redden- 
ing. However, there is evidence, mainly from the less strongly 
polarized objects, that a young stellar population can make a 
significant contribution to the UV continuum. 
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Object 


z 


P(%) 


P 


Lya/ClV 


Nv/ClV 


Hen/Civ 


FSC 10214+4724 


2.282 


26±2 




0.276±0.003 


0.736±0.001 


0.496±0.001 


4C+41.17 


3.798 


<2.4 (2a) 


-1.8 


11.1±0.7 


0.29±0.03 


0.42±0.02 


3C 256 


1.824 


10.9±0.9 


-0.71 


10.4±0.2 


0.27±0.17 


1.05±0.02 



Table 4. Polarization and line ratio results from the literature for three sources which have similar quality spectropolarimetry. 
The data are taken from: FSC 10214+4724, Goodrich et al. (1996); 4C+41.17, Dey et al. (1997); 3C 256, Dey et al. (1996), 
Simpson et al. (1999). The fractional polarization measurement has been chosen to match as closely as possible the aperture and 
waveband we have used for our sample sources. In the case of 3C 256, no polarization measurement has been made shortward of 
Civ, but Dey et al. (1996) show that P is quite flat from 1500-2500 A. When forming line ratios, we have excluded components 
with a FWHM>3,000 km s" 1 . 



4.1.1. The scattered continuum 

The high continuum polarization shows that a large fraction of 
the UV continuum must be due to scattered AGN light in all 
objects in our sample except possibly for the least polarized 
source (0731+438). This case is strengthened by the clear de- 
tection of broad lines in the most polarized objects (021 1 — 122, 
4C+23.56) with equivalent widths which are similar to those 
seen in quasars. 

The remarkable similarity of the shape of the HzRG ultra- 
violet spectral energy distribution (SED) with that of quasars 
shows that the scattering process must be approximately inde- 
pendent of wavelength. The most natural explanation for such 
'grey' reprocessing is that the dominant mechanism is Thomp- 
son scattering by electrons. However, this interpretation en- 
counters several difficulties. Scattering by a population of hot 
(T e > 10 7 A") electrons is ruled out because the Doppler broad- 
ening would smear out the scattered broad lines. In the two ob- 
jects in our sample for which we can easily measure the width 
of the scattered broad CIVA1550 line we can estimate an up- 
per limit on T e of ~ 10 5 A (see Miller, Goodrich and Math- 
ews 1991 and Cimatti et al. 1996) assuming a typical value of 
~ 6500 km s -1 for the width of the CIVA1550 line of a ra- 
dio loud quasar (as measured on Cristiani and Vio 1990 radio 
loud quasar composite spectrum). In addition, estimates of the 
expected X-ray fluxes for such a population of hot electrons 
in distant 3CR radio galaxies are in strong disagreement with 
observations (see for instance calculations made by Dey et al. 
1996 for 3C 256). Scattering by a cooler (T e ~ 10 4 - 10 5 K) 
population of electrons as the dominant mechanism is more dif- 
ficult to rule out. The most compelling argument is based on the 
far greater scattering efficiency per unit mass of dust than elec- 
trons which implies that, as soon as some dust is present in the 
ISM, it will dominate the scattering process, an assumption we 
maintain throughout the rest of this discussion. 

Standard homogeneous optically thin dust scattering mod- 
els predict a bluening of the emergent spectrum due to the 
wavelength dependence of the scattering efficiency of small 
dust grains. In order to reproduce observations, previous mod- 
eling has been based on optically thin scattering followed by 
arbitrary extinction by a screen of dust. The implicit assump- 
tion is that some part of the dusty medium is illuminated by the 
quasar and scatters light to the observer while another part is 
in the 'shadow' and produces extinction of the scattered light 



(see Manzini and di Serego Alighieri 1996). A more realistic 
model employs the same population of dust grains to scatter 
and absorb the light whereupon the emergent scattered light 
is maximized where a photon experiences scattering and ex- 
tinction optical depths close to unity (j scat ~ r ext ~ 1 ). If 
the medium is inhomogeneous (clumpy), a natural luminosity 
weighting process operates by which the emergent flux over a 
wide range of wavelength comes from paths with optical depth 
close to one. The scattered spectrum is then approximately the 
input spectrum multiplied by the dust albedo (the ratio of the 
scattering to the extinction efficiency) which, for Galactic dust 
is almost grey at ultraviolet and optical wavelengths. 

Radiative transfer in a clumpy dusty medium has been ex- 
tensively studied (see eg. Witt and Gordon 1996,1999; Varosi 
and Dwek 1999 and references therein) and applied to sev- 
eral types of astrophysical objects such as circumstellar shells 
and starforming galaxies. The expected greyness of the scat- 
tering over a large wavelength range described above is con- 
firmed by all studies. We compared our observations to re- 
sults obtained using the analytical model developed and teste d 
against Monte-Carlo simulations by Varosi and Dwek ( 1999 ). 
These analytical approximations, based on an extension of the 
mega-grain model of a two -phase medium first proposed by 
Hobson and Padman ( 1993 ) allowed us to easily explore the 
main parameters describing the clumpy medium: the clump 
volume filling factor f c and the ratio of the cloud to inter- 
cloud medium density a (see Varosi and Dwek 1999 for de- 
tails). In our simulations, we used a 40 kpc diameter spheri- 
cal two-phase dusty medium illuminated by a q uasar central 
source with spectrum from the Cristiani and Vio (1990) quasar 
composite. The adopted total dust mass is 10 8 M Q typical of 
dust masses inferred from submillimeter observations (Best et 
al. 1998a, Cimatti et al. 1998a, Archibald et al. 2000) and we 
used the standard Galactic dust model by Mathis, Rumpl and 
Nordsiek (1977). The result is displayed in fig. || where the 
HzRG composite spectrum (grey line) is compared to the scat- 
tered quasar (black continuous line). This provides a remark- 
ably good fit to both continuum and broad lines. We also plotted 
the original quasar composite (dotted black line) as a reference 
to illustrate the approximate greyness of the scattering. 

The best results were obtained for models with high clump 
to inter-clump medium density contrast (a — 1000) and low 
clump volume filling factor (/ c = 0.01). The reason for this 
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Fig. 5. Comparison of the HzRG composite spectrum (in grey) with quasar radiation scattered by a clumpy medium (continuous 
black line) modeled using Varosi and Dwek (1999) models with high density contrast (a — 1000) and small clump volume 
filling factor f c — 0.01. This medium produces an almost grey scattering illustrated by the similarity between the scattered 
spectrum and the input quasar composite spectrum (dotted line). The difference between the quasar and the scattered spectrum 
around 2200 A is due to the dust albedo feature. Note how well the HzRG composite spectrum tracks this feature. Spectra are 
normalized at 2000 A. 



is that in order to have grey scattering, high density clumps in 
an extended low density inter-clump medium are required. In 
this regime radiative transfer is dominated by the dense clumps. 
Additionally, this ensures that the direct view to the quasar is 
usually relatively unobscured, consistent with the scarcity of 
highly reddened radio loud quasars (see eg, Simpson and Rawl- 
ings 2000). However it is interesting to note that this scenario 
also predicts that there must be cases where our line of sight to 
the quasar intercepts one of these high density clumps, in which 
case the direct view to the quasar would be completely blocked. 
Scattered light would then dominate the flux in the UV-optical, 
resulting in a type 1 object masquerading as a radio galaxy. 

The main features of the continuum shape — the dip around 
2200 A and the rollover below 1500 A — are well reproduced 
by this model (see fig. ||). It is difficult, however, to use these 
observations to reliably constrain the dust properties since 
the scattered spectrum is sensitive only to the difference in 
wavelength dependence of the scattering and absorption cross- 
sections. In particular the possible detection of the 2200 A dust 
feature is difficult to confirm because it overlaps with the UV 



Fell blends often present in quasar spectra above 2300 A (the 
so-called small blue bump). In addition, the signal-to-noise ra- 
tio in this region of the spectrum is rather low due to strong sky 
lines (shaded areas on fig. [[]). The strength of these Fells blends 
in quasar spectra varies significantly from object to object and 
its origin is still poorly understood. It seems, nevertheless, diffi- 
cult to ascribe the observed feature to Fell blends alone since it 
would require this blend to be unusually strong in the spectrum 
of the hosted quasar in all sources, as strong as that observed 
in ultra-strong Fell emitters like quasar 2226—3905 (Graham 
et al. 1996). In consequence, the feature we observe around 
2200 A is probably a combination of the 2200 A dust feature 
and UV Fen blends above 2300 A. 

Although inhomogeneous radiative transfer models consid- 
ering anisotropic central sources are not yet available, we can 
check that the predicted scattered fluxes are consistent with our 
observations. The model described above predicts that the frac- 
tion of scattered light in a spherical medium is of the of order 
of 7% of the input flux from the central source between 1000 
and 3000 A. This value is consistent with the Witt and Gor- 
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don (1999) "clumpy shell" model predictions (their t v = 0.75 
model). Assuming an average half opening angle of 45° for 
the obscuring torus, we then estimate that about 2% of the cen- 
tral AGN flux is scattered by the ISM of a radio galaxy. The R 
magnitudes of the objects in our sample range between 21 and 
23 and the inferred magnitude of the central source is there- 
fore between 17 and 19, consistent with observed bright quasar 
magnitudes at z ~ 2.5. 

The spatially-integrated level of continuum polarization 
is determined by two processes: the geometrical dilution — 
which is the result of averaging over the possible scattering an- 
gles within the ionization cones — and dilution by sources of 
unpolarized radiation. The appropriate biconical geometry in 
the framework of the unified scheme has be en studied in detail 
by Manzini and di Serego Alighieri ( 1996| ) assuming an opti- 
cally thin homogeneous medium. Currently, none of the avail- 
able radiative transfer models for inhomogeneous clumpy me- 
dia have considered the fractional polarization of the scattered 
radiation in a biconical geometry. However, a study of opti- 
cal depth and ge ometrical effects on dust scattering by Zubko 
and Laor ( 1999 ) shows that the fractional polarization is al- 
most independent of optical depth due to the relatively small 
contribution of multiple scattering to the net polarization. In 
consequence, we assume that the calculations made for ho- 
mogeneous optically thin medium by Manzini and di Serego 
Alighieri (1996) provide a reasonable estimate of the fractional 
polarization of the scattered radiation. Considering only calcu- 
lations made for standard MRN interstellar dust composition 
and a torus half opening angle $ = 45°, the fractional polar- 
ization at a given wavelength is determined by the observer's 
viewing angle 9. These models predict that the continuum po- 
larization at 1500 A varies between 25% for 9 = 90° (when the 
torus is seen edge on) and about 10% for 9 = 50°. This range in 
fractional polarization depends on the assumed torus opening 
angle. The value we use in this calculation is illustrative and is 
consistent with Barthel's (1989) statistical analysis. 

The range in observed continuum polarization in our sam- 
ple is compatible with these predictions. It goes from less than 
3% in 0731+438 to about 20% in 0211-122 longward Lya, 
which is lower than the range spanned by models. Using the 
following equation describing the dilution of polarization of 
the scattered radiation F scat of fractional polarization P sca t by 
an unpolarized component F unpo i : 



100 



Pobs — P. st 



Fscat + F un p i 



we can set limits to the contribution of unpolarized radiation 
at ^1500 A making extreme assumptions for the unknown ob- 
server's viewing angle as displayed on fig. ^(see also table |]). 
While the polarized scattered radiation must contribute more 
than 80% in 02 1 1 — 1 22, unpolarized radiation constitutes more 
than 70% of the continuum in the least polarized radio galaxy in 
this sample. All objects except 0731+438 are compatible with 
scattered continuum accounting for at least half of the observed 
continuum. The two main sources of unpolarized radiation — 
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Fig. 6. Constraints on the fraction of unpolarized light con- 
tributing to the continuum at -1500 A. The two lines delin- 
eating the shaded area represent F unpo i for extreme values of 
the observer's viewing angle (top line: 9 = 90°; bottom line: 
9 = 50°). Vertical labeled lines represent the possible range of 
variation of F unpo i for each object in our sample. 



the nebular continuum and a young stellar population 
discussed in the next two sections. 



4.1.2. The nebular continuum 



As shown by Dickson et al. (1995), the nebular continuum 
makes a significant contribution to the UV continuum of pow- 
erful radio galaxies. Its spectrum, normalised to the observed 
Hj3 flux, is computed using tabula ted coefficients for T e = 
10 A K given in Aller et al. ( |l987| ). Since the Lya flux is 
strongly influenced by neutral hydrogen as well as dust, the 
Lya/H(3 ratio does not provide a reliable estimate of the Hp 
flux. We use instead the Hell line flux and th e average value 
Hell/H/3 ~ 3.18 given by McCarthy ( |l993| ) to compute the 
contribution of the nebular continuum. Assuming that both the 
emission lines and the nebular continuum are affected by the 
same amount of dust extinction, our estimate based on the Hell 
flux approximately takes this into account. The contribution to 
the total flux at 1500 A for each object is given in the third col- 
umn of table U In general, the nebular continuum contributes 
less than 25% of the continuum at 1500 A leaving room for a 
significant starlight contribution of unpolarized radiation. 

4.1.3. The young stellar population 

From the constraints on scattered light and estimates of the 
nebular continuum, we obtain limits on the fraction of contin- 




Fig. 7. Three components fit to the two objects in our sample showing extreme polarization properties. Left: 0211 — 122. Right: 
0731+438. 



Object 


Scattered 


Young stellar 


Nebular 




AGN (%) 


population (%) 


continuum (% 


0211-122 


80-91 


11-0 


9 


4C+23.56 


61-93 


32-0 


7 


4C-00.54 


47-88 


41-0 


12 


4C+03.24 


45-97 


52-0 


3 


0828+193 


40-75 


35-0 


25 


4C+48.48 


30-74 


58-14 


12 


4C+40.36 


29-73 


61-17 


10 


0943-242 


27-66 


53-14 


20 


0731+438 


0-30 


84-54 


16 



Table 5. Range in the relative contribution of the different con- 
tinuum component at 1500 A. 



uum contributed by a young stellar population. These limits 
are shown in the second column of table |[ In all objects except 
0731+438, starlight contributes less than half of the UV con- 
tinuum. In consequence, it is not possible to accurately charac- 
terise the underlying young stellar population. 

Unambiguous direct spectral signatures of young massive 
stars that we might expect to detect are SilV A 1400, Civ A 1549 
P Cygni profiles and purely photospheric absorption lines (ie. 
uncontaminated by interstellar absorption) from O and B stars 
like Sim A 1294, Cm A 1427, Sv A 1502. We have searched for 
the presence of such features in our sample, especially in low 
polarization objects where the contribution of a young stellar 
population is likely to be the most important. 

We find no clear evidence for the presence of stellar 
P Cygni profiles in any of our spectra. These features are 
prominent in very young stellar populations (t < 5 Myr) 
with equivalent width W\ ~ 8 A (as measured on Civ line 
in Leitherer et al. 1999 [Starburst 99] 1 Myr solar metallic- 
ity Salpeter initial mass function (IMF) instantaneous starburst 
model). Such features are prominent enough to be detected 
given the quality of our data. However, the equivalent widths 
of these features vary strongly with metallicity by a factor of 



about 5, the more metal-rich starbursts having stronger lines 
(Heckman et al. 1998). In addition, the presence of strong nar- 
row emission lines and scattered broad lines could, in some 
objects, partially fill such features if present. The blue wing 
of the Civ A 1549 line profile is also affected by the strong 
Sill AA 1527, 1534 absorption and emission line. 

Stellar population synthesis models and observations of 
both local starbursts (eg. Leitherer et al. 1996 ; Conti et 
al. 1996) and high redshift star forming galaxies (eg. Pettini 
et al. 2000) show that the equivalent width of Sim A 1294, 
Cm A 1427 and Sv A 1502 photospheric absorption lines varies 
from 0.6 to 0.8 A. The strength of these lines also depends 
strongly on metallicity. Heckman et al. (1998) can only ten- 
tatively detect photospheric features in their high signal-to- 
noise low metallicity composite spectrum. We estimate the 
detection limit for an unresolved absorption line to be about 
W\ ~ 0.7 A in the restframe, precisely in the expected range 
of observed values for stellar features. This means that in our 
observations we could only hope to detect these absorption 
lines in the most favorable cases where young stars dominate 
the continuum. Non-detections do not, therefore, provide use- 
ful constraints. We do detect an absorption line that could be 
identified as Cm A 1427 photospheric line in our spectrum of 
0731+438 (and also in 4C— 00.54, see discussion in paper I). 
The measured equivalent width in the restframe is 0.85+0.2 A, 
consistent with values observed in starbursts. This purely pho- 
tospheric line is present only in the spectrum of early-B stars 
(de Mello et al., 200C ) but since no other characteristic feature 



of young massive stars is detected in any of our spectra, we 
consider this identification as very uncertain. 

Constraints on the star formation rate (SFR) can be derived 
by comparing the luminosity at 1500 A (listed in table |2|), cor- 
rected for scattered light and nebular continuum contributions 
(L* 5Q0 ), to predictions from stellar population synthesis mod- 
els. We used a Starburst 99 (Leitherer et al. 1999), solar metal- 
licity, Salpeter IMF, continuous star formation model. In our 
sample, log(L* 500 ) ranges from about 42.2 to less than 40.7. 
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Assuming that the stellar population is older than 10 Myr (ie. 
when the stellar population has reached an equilibrium after 
the initial onset of the burst), this translates into a range in SFR 
from -60 Mqi/t^ 1 to less than 2 M^yr^ 1 , not taking into ac- 
count any dust reddening correction. 

We could, in principle, combine the luminosity informa- 
tion with the spectrum colour to constrain both the redden- 
ing and the age of the stellar population. However, since re- 
processed AGN radiation dominates the continuum of all ob- 
jects in this sample with the exception of 0731+438, the slope 

does not provide meaningful information about the proper- 
ties of the stellar population. In this one exception, the slope at 
1500 A ~ —1.4, is significantly redder than values spanned 
by starbursts models for which ranges from -2.5 at 10 Myr 
to -2.0 at 1 Gyr. This shows that in this object the stellar popu- 
lation is reddened with Eb-v values between ^0.1 and ^0.2 
if the age of the stellar population (t*) is between 10 Myr and 

1 Gyr. The range in L^qq for this object then translates into 
a SFR between 30 and 120 M^r- 1 for E B -v = 0.1 and 
t* = 10 V or between 60 and 260 Mq?/^ 1 for E B -v = 0.2 
and t* — 10 7 yr using a Galactic reddening law (consistent 
with our scattering model) which gives a dust extinction in 
magnitudes of A1500 = 8.3 x Eb-v- 

This range in SFR appears to be comparable to what is 
typically observed in Lyman break galaxies at a similar red- 
shift. For instance Sawicki and Yee 1998 find that the me- 
dian dust-corrected SFR in z < 3 Lyman break galaxies is 
~ 560Af Q i/r~ 1 (with a typical reddening Eb-v ~ 0.3) or 
~ 20M & yr~ 1 without dust correction. In the well studied 
lensed z ~ 2.7 Lyman break galaxy MS 1512— cB58, Pettini et 
al. (2000) derive a SFR of about 70A/ Q yr _1 (all values were 
converted to the cosmology used in this paper). 

If we now compare to instantaneous burst models, we can 
set limits on the total mass A/* and the age t* of the burst 
in 0731+438, depending on reddening. We obtain 4.9 10 9 < 
M* < 2.210 10 M o and t* ~ 4.7 10 7 yr without redden- 
ing, 2.6 10 9 < M* < 1.210 10 Af Q and t* ~ 1.6 10V 
for E B -v = 0.1 and 5.6 10 8 < M* < 2.51O 9 M and 
t* ~ 0.12 10V for E B -v = 0.2. 

As an example we present in fig. ^ a three component fit 
(scattered quasar, nebular continuum and young stellar popu- 
lation) of the two objects in our sample showing extreme po- 
larization properties. The contribution of each component for 
these objects is given in table |[ For 0211 — 122, — the most 
polarized object — we display the fit that maximizes the stel- 
lar contribution allowed by our model. On the other hand for 
0731+438 — the least polarized object — we show the so- 
lution which maximizes the scattered quasar's contribution. In 
both cases the stellar population was modeled using the Star- 
burst 99 solar metallicity 20 Myr instantaneous burst model. 

The conclusion of this analysis is that the brightness of the 
scattered quasar light in these powerful HzRG can swamp a 
SFR of several tens of Solar masses per year unless a burst is 
observed during its initial, very luminous phase. One at most 
of the objects in our sample appears currently to be undergoing 
such a burst. 



4.2. Predicted far infrared continuum 

Current submillimeter instruments are just reaching the sensi- 
tivity required to start measuring the far infrared (FIR) emis- 
sion of brightest HzRG and provide us with useful flux upper 
limits (Archibald et al. 2000). As a consistency check, we have 
computed the FIR SED that is expected from HzRG in a way 
that is consistent with our dust scattering model. This estimate 
is restricted to the FIR flux re-radiated by the same dust which 
is responsible for that scattering of the ultraviolet continuum 
(ie. within the torus opening angle). Since much of this dust is 
located at large radial distances from the nucleus, it is relatively 
cool. 

The characteristics of the ISM are the same as that used 
to model the scattered continuum: a 40kpc sphere containing 
10 8 A/q of Galactic type dust in an clumpy geometry with a 
density contrast a = 1000 and filling factor f c — 0.01 is illu- 
minated by an anisotropic cent ral sou rce. The method is fully 
described in Varosi and Dwek ( 1999). 



This model does not provide a full treatment of the effect 
of the anisoptropy of the central source nor does it include the 
effect of radiation below the Lyman limit (-30% total energy 
in the UV). The latter would require modeling the ionization of 
gas and the heating of dust by absorption of Lya photons that 
is beyond the scope of this paper. 

Under these assumptions, the result of the simulations is the 
following energy balance. The total input energy is ~ 10 14 Lq 
— the integrated quasar ultraviolet-optical continuum above 
the Lyman limit. We are only interested here in the fraction of 
this energy (29%) that escapes within the torus opening angle 
(the remaining 7 1 % is absorbed and re-radiated by the dusty 
torus at wavelength around 10 /Ltm). The ^3 10 13 L Q going 
into the illuminated cone then divides into 77% escaping di- 
rectly without interacting with dust grains, 7% that is scattered 
(see sect. |4.1.l| ) and 16% absorbed and re-radiated in the FIR. 
The resulting FIR luminosity of the dust within the opening an- 
gle of the torus is about 5 lO 12 ^©. More than 75% of this en- 
ergy comes from the clumps and the rest from the interclump 
medium. 

In order to compute the FIR SED, the temperature equilib- 
rium for graphite and silicate grains is computed by equating 
absorbed and emitted luminosities in clumps and in the inter- 
clump medium using a power law temperature probability dis- 
tribution and assuming that that the radiation density decreases 
with distance from the source like an inverse powerlaw with 
index 2.5 (2 would be for optically thin case, see Varosi and 
Dwek 1999). 

The resulting modeled SED from 0.1 to 1000 /j,m is plot- 
ted in fig. ||. The calculation shown here has been made to 
match the scattered flux of 4C+48.48 from our Keck spec- 
trum because this is the only object in our sample detected with 
SCUBA at 850 fim. This corresponds to a quasar of magnitude 
R = 17.3. The scattering is almost grey in the UV but at longer 
wavelengths, clumps become optically thin and the scattering 
deviates from grey and eventually the medium becomes trans- 
parent (the radiation that escapes directly [green curve in fig. 
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Fig. 8. HzRG SED from 0. 1 to 1000 [im. Black continuous line: predicted scattered (left) and FIR (right) continuum; black dotted 
line: input UV-optical quasar spectrum; green line: radiation escaping directly; red line: composite HzRG spectrum scaled to 
4C+48.48 luminosity. Dots and arrows show SCUBA measurements and 2<r upper limits at 450 [im and 850 fim (light blue) of 
HzRG between redshift 2 and 3 from Archibald et al. 2000; red dot and arrow: SCUBA measurements for 4C+48.48 at 850/im 
and 450/jto. 



|8p converges toward the input spectrum [dotted curve] and the 
scattered radiation drops longward of 1/im). 

Our predictions are approximately a factor of ten below the 
upper limits. It is important to stress at this point that what 
we compute here is clearly a lower limit to the FIR flux for 
several reasons. The contribution of the dusty torus IR radiation 
is not included but it will mainly contribute around 10 (im. In 
this calculation only the radiation of the dust within the cone 
opening angle is computed. The rest of the host galaxy (-70% 
of the ISM) mainly heated by star formation and reprocessed 
AGN radiation is not taken into account although it clearly has 
a significant contribution at long wavelength (A > 100 /Ltm). 
Also, as mentioned above, the transfer of the radiation below 
the Lyman limit is not taken into account. 

The point of this calculation is not to obtain the best fit to 
the data but to show that our model, which is consistent with 
the observed scattered UV continuum, does not over-predict 
the FIR radiation. This confirms the impression we have from 
observations that current instrumentation is really at the sensi- 
tivity limit needed to detect HzRGs and that we currently just 
pick the brightest objects. Most powerful radio galaxies at red- 



shift between 2 and 3 should be detected in the near future with 
an increase of instrumental sensitivity less than a factor of ten. 

4.2.1. Emission lines 

The emission lines seen in radio galaxies arise from two quite 
distinct regions. We see scattered light from the central quasar 
which is characterised by its polarization and includes both 
broad and narrow components. These narrow components arise 
in a region which is larger than the broad line region (BLR) 
proper, but still small enough to be obscured from our direct 
view, with the lines coming from transitions with a high criti- 
cal density for collisional de-excitation being more prominent 
(Hes, Barthel & Fosbury 1993). The emission lines from the 
global ISM of the host galaxy come predominantly from those 
zones which see the ionizing radiation from the quasar directly 
and are generally much stronger that the scattered narrow lines. 
While a number of mechanisms may be responsible for heat- 
ing and ionizing this gas, including shocks and cosmic rays, 
the ionization of the warm (~ 10 4 K) ISM in the hosts of these 
powerful radio galaxies is thought to be dominated by the AGN 
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radiation field (Villar-Martfn et al. 1999). The presence of this 
galaxy-scale gas, illuminated by an intense ionizing radiation 
field resulting in bright emission lines spanning a broad range 
of ionization states, presents an opportunity for detailed stud- 
ies of its physical state and chemical composition which is 
absent in galaxies without AGN. This opportunity cannot be 
easily exploited in the quasars themselves since their spectra 
are so dominated by direct light from the AGN. It is in the ra- 
dio galaxies that we can most effectively exploit the 'natural 
coronograph' provided by the local obscuration of the quasar. 

Studies of chemical composition are likely to be particu- 
larly revealing since this is the epoch where we expect the mas- 
sive host galaxies and the AGN themselves to be in the process 
of assembly or at least to be in an early stage of their evolution. 
Since the timescales for AGN phenomena and galactic evolu- 
tion are very different — by factors of a hundred or so — we 
can consider the quasar as acting like a 'flash bulb' which gives 
us a snapshot view of its host galaxy. By studying a number of 
such objects at a similar redshift, we might expect these snap- 
shots to reveal hosts in different stages of evolution even if the 
galaxy and AGN formation triggers are closely related. 

In table ^, we present measurements of all the detected 
emission lines. The spectrograph aperture ranges from 4 — 8" x 
1" which corresponds to 500-1000 kpc 2 at these redshifts. For 
objects of this size and complexity, we must be realistic in our 
expectations for the results of the analysis. The volume inte- 
grals represented by emission line measurements are taken over 
a large fraction of the ISM of a whole galaxy and will include 
inhomogeneities on many scales. The modelling we do vastly 
oversimplifies the situation but nonetheless may allow us to in- 
terpret differences in line ratios from object to object in terms 
of variations in their properties even if we can place less confi- 
dence in, for example, absolute elemental abundance determi- 
nations. 

The spectra are characterised by a high state of ionization 
with the Civ, Nv and, in 0943-242 which has a sufficiently 
high redshift, Ovi resonance lines being prominent. The Hell 
recombination line is also strong, comparable with Civ. The 
signal-to-noise of the spectra is sufficient to measure a number 
of intercombination lines from carbon, nitrogen, oxygen and 
silicon. The average relative intensity of the strong lines is sim- 
ilar to that in the average MRC radio galaxy spectrum compiled 
by McCarthy (1993). The most significant variation amongst 
the stronger lines is seen in Lya, which is clearly affected by 
self-absorption, and in Nv. 

4.2.2. Photoionization modelling 

A set of diagrams containing the line ratio data from our sam- 
ple is shown in fig. ^ and [fo| The behaviour of the first of these, 
the Nv/Hell vs. Nv/Civ diagram (fig. ||), is reminiscent of the 
behaviour of luminous QSO broad emission line regions dis- 
cussed by Hamann and Ferland ( |l993| ; |l999| ) (HF93, HF99). 
These authors showed that high redshift quasars (z > 2) define 
a tight correlation in this diagnostic diagram and the model- 
ing of the emission line ratios led them to conclude that the 



two Nv ratios imply supersolar metallicities in the broad line 
region of many high redshift quasars. They interpret the cor- 
relation in the Nv diagram as a sequence in metallicity such 
that the highest redshift/most luminous objects show the high- 
est metallicities (up to lOx Zq or so). The reason for this could 
be related to the higher QSO and/or host galaxy masses at large 
redshifts, implying a mass-metallicity relationship in QSO sim- 
ilar to that seen in nearby ellipticals (Tinsley 1980; Pagel & 
Edmunds 1981; Vader 1986; Bica 1988 ). This diagram is fur- 
ther interepreted in terms of elliptical galaxy chemical evolu- 
tion models by Matteucci & Padovani (1993). 

The correlation defined by the sample of HzRG suggests a 
similar interpretation: we might be witnessing different levels 
of metal enrichment of the gas from object to object extend- 
ing, possibly, to supersolar metallicities. If indeed we are see- 
ing large metallicity variations within our sample, the fact that 
it is apparent in the narrow emission lines emitted by the galaxy 
ISM provides a much stronger argument for global galactic 
chemical evolution than the small quantity of nuclear gas rep- 
resented by the quasar BLR. 

In order to test the validity of this interpretation it was first 
necessary to explore whether models other than a metallicity 
sequence could reproduce the observations. Villar-Martfn et al. 
(1999) studied the effects of shock ionization, AGN photoion- 
ization and the influence of the AGN continuum shape, density 
and/or ionization parameter They showed that these models 
could explain neither the Nv correlation nor the very strong Nv 
emission observed in some objects. Other effects investigated 
are the depletion of carbon onto dust grains and the effect of 
illumination/viewing geometry on resonance line transfer (see 
Villar-Martfn, Binette & Fosbury 1996 ). Carbon depletion can 
produce large Nv/Civ ratios, but Nv/Hell remains nearly con- 
stant unless the metallicy is varied. Regarding geometrical ef- 
fects, since Nv and Civ are resonance lines, the contribution 
relative to other emission lines should vary with the orientation 
of the object. However, Nv/Civ and Nv/Hell vary in opposite 
ways with perspective (Nv/Civ increases when Nv/Hell de- 
creases and viceversa) so that it is not possible to produce large 
values (>1) of Nv/Hell and Nv/Civ at the same time. 

The effect on the Nv behaviour of a metallicity sequence 
was then examined to see if it could explain the observations 
in this diagram and remain consistent with the other emission 
line ratios. Thanks to the high signal-to-noise of the spectra, 
we could use many fainter emission lines to to compare with 
the models. It was assumed that the gas (100 cm~ 3 ) is pho- 
toionized by a power law of index a — —1.0 (Villar-Martfn et 
al. 1999) and the same U(= 0.035) was used for all the ob- 
jects. Such a constant ionization state is suggested by the small 
variation of Civ/Cm] and Civ/Hell, fig. |To| ). 

It is found that: 

- A sequence in metallicity can reproduce both the observed 
correlation and the strength of the Nv emission. The heavy 
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element abundances relative to H vary between 0.4 and 
4-Zq in the models. 

- The N abundance must increase quadratically instead of 
linearly with respect to carbon, oxygen etc. 

- There is good agreement between the model predictions 
and the data in all diagrams except for Niv]A1488. This 
line is predicted to be stronger than observed. None of 
the models we explored can explain this discrepancy. A 
similar inconsistency was reported for the Seyfert galaxy 
NGC1068 (Kraemer & Crenshaw 2000). However, the fact 
that both the data and the models define a tight correlation 
in the Niv]A1488 diagram (see also Olll] diagram) sup- 
ports a metallicity sequence. 

Another possibility, based on the polarization results, is that 
the extended gas is mixed with dust. We studied the effects 
of internal dust as a variation of the dust to gas ratio (deple- 
tion, scattering, absorption) and found that not only are high 
metallicities needed, but also unrealistically high densities (10 6 
cm -3 ) in the extended gas to explain the relative strength of N V 
in some objects. In addition, these models produce strong dis- 
crepancies with other observed UV line ratios. These results 
imply that the emission lines come preferentially from regions 
where the effects of internal dust are small. 

We conclude that the Nv diagram implies solar or super- 
solar metallicities in the extended gas of many HzRG. The se- 
quence shows that the level of enrichment varies from object to 
object. The quadratic increase of N abundance suggests that the 
N production is dominated by secondary nucleosynthesis pro- 
cesses, as expected for high metallicities (Henry et al. 2000). 

4.2.3. Correlations 

Amongst the data for the continuum and the stronger emission 
lines, there are two correlations which are immediately appar- 
ent. The first is between the continuum fractional polarization 
P and the strength of Lya. Fig. [IT] shows that the highly po- 
larized sources have low Lya/Civ ratios and fig. [j], especially 
the spectrum of 0211-122, suggests that the weakness of Lya 
is due to self absorption as has b een de monstrated to occur in 
radio galaxies by van Ojik et al. (1997a). 



The second correlation, shown in fig. |l2|, is between P 
and the Nv/Civ ratio with the latter spanning a factor of five 
within our sample. The fact that this ratio is uncorrected with 
Clll]/Civ or Hell/Civ suggests that this is not simply an ion- 
ization effect. The correlation is also present between Nv/Hell 
and P, which suggests — as discussed in the previous section 
— that carbon depletion onto grains is not the dominant effect. 
The tightness of this correlation is remarkable since, whatever 
the underlying cause, these two physical quantities are likely to 
represent very different averages over the emitting volumes. In 
the simplest case, the line ratio is a microphysical property of 
the ionized gas while the continuum polarization is dependent 
on large-scale structure and orientation. 
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Fig. 11. The relationship between the strength of Lya — ex- 
pressed as a ratio to Civ — and the continuum fractional po- 
larization measured in the bin just longward of Lya. The error 
bars represent ler statistical uncertainties in P but, for the line 
ratios, are derived from uncertainties in the continuum fitting 
(see text). For values of P < 3%, the points are plotted as 3ct 
upper limits. 
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Fig. 12. The relationship between the narrow line ratio Nv/Civ 
and the continuum fractional polarization measured in the bin 
just longward of Lya. 
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Fig. 9. The Nv/Hell vs. Nv/Civ diagram as plotted and modelled by HF93 for the BELR of QSOs — which clearly have much 
higher densities than the gas we are sampling in the narrow lines. The dashed line represents the locus of one of their chemical 
evolution models: M4a which covers a range of 2-10 times solar for the primary elements in this frame. The sequences implies 
solar or supersolar metallicities in the extended gas of many HzRG and different levels of enrichment in different objects. The 
three objects from the literature (see table 0) are plotted with different symbols (triangles). 
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Fig. 10. Diagnostic diagrams involving some important UV emission line ratios including the intercombination lines Niv] A 
1486, Om] AA 1661,1666 and Cm] A 1909. These include only data points from our sample. The model metallicity sequence 
that best reproduces the data trend (solid line) shows good agreement with the data in most diagnostic diagrams although the 
Niv] line predictions are substantially displaced (lower left, see text). The numerical labels show the metallicity sequence from 
0.4 to 4 times solar for elements other than nitrogen which increases quadratically. The top two frames show the small variation 
in ionization state within the sample while the oxygen diagram (lower right) confirms the existence of the metallicity sequence. 
In this and the associated figures, only emission line components with a FWHM of < 3, 000 km s _1 are included. 
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The nature of these correlations can be further illustrated by 
comparing our sample with other objects in the literature which 
have comparable quality spectroscopy (see table Q). Our most 
polarized source, 0211 — 122, has a line and continuum spec- 
trum which is remarkably similar to the lensed ultraluminous 
IRAS galaxy FSC 10214+4724 (z = 2.282, Goodrich et al. 
1996). Both sources have weak (absorbed) Lya, strong Nv and 
a high continuum polarization: 20% for 0211 — 122 and 26% 
for FSC 10214+4724 measured in the same restframe band 
between Nv and SilV. Although FSC 10214+4724 is not a 
powerful radio source, the Keck spectropolarimetry shows it 
to contain a luminous hidden quasar. 

Our lowest polarization source, 0731+438, has strong Lya 
and weak Nv similar to 4C+41.17 (z = 3.798) which Dey et 



al. (1997) have argued has a continuum dominated by a young 
stellar population. The relatively strong, low ionization inter- 
stellar absorption lines seen in the 4C object appear, however, 
in emission in 0731+438. 

The z — 1.824 radio galaxy 3C 256 has both imaging- 
(Jannuzi et al. 1995 ) and spectropolarimetry (Dey et al. 1996) 
with P ranging from 1 1-18% depending on aperture and wave- 



length. Palomar 5m spectroscopy by Simpson et al. ( 1 999 ) pro 



vides measurements of the Lya, Nv, Hell and Civ emission 
line fluxes. This object has a continuum shape between 1500- 
2600 A which is essentially identical to our average spectrum 
and is well fitted by the reflected quasar model. 

These three sources have been plotted with symbols differ- 
ent from the main sample in the figures. 

In order to understand the origin of these correlations we 
need to review the possible physical reasons for the variations 
of P, Lya and Nv. For the scattering of quasar light by dust, the 
fractional polarization we measure is determined by the partic- 
ular biconical scattering geometry and global orientation and 
by the presence of diluting, unpolarized continua. The strength 
of Lya with respect to the other emission lines is probably most 
strongly influenced by neutral hydrogen absorption with the 
possibility of dust playing a role in quenching multiply scat- 
tered line photons (Villar-Martm, Binette and Fosbury 1996). 

We have concluded that the variation of the Nv/Civ ratio 
is the result of abundance variations. We cannot currently pro- 
vide a detailed explanation for this behaviour but it is likely, 
we believe, to be a manifestation of an evolutionary sequence 
where a host galaxy starts its assembly with low metallicity 
and being relatively dust-free. As the spheroid building pro- 
ceeds, the chemical enrichment of the ISM is accompanied 
by dust production and dispersal. The similarity between the 
spectra of 0211 — 122, our most polarized HzRG, and several 
of the AGN-powered ULIRG, eg. FSC 10214+4724, suggests 
that these dust-enshrouded objects are the end point of a Gyr of 
spheroid assembly. 



5. Conclusions 

By performing spectropolarimetry of the restframe ultraviolet 
emission from a sample of z ~ 2.5 radio galaxies, we have 



been able to isolate and identify the dominant contributors to 
their observed spectra. These are: 

- Spatially-extended, dust-scattered light from the underly- 
ing obscured quasar. Due to the supposed dumpiness of the 
scattering medium, a natural 'luminosity-weighting' pro- 
cess operates which means that the scattered photons orig- 
inate predominantly from regions where r ~ 1. This re- 
sults in approximately grey reflection of a few percent of 
the quasar luminosity. Deviations from greyness around 
2200 A are expected for standard Galactic dust and are 
seen in our data if the quasars do not have abnormally 
strong Fell emission. Based on the spectral shape and the 
fractional polarization, the scattered light accounts for be- 
tween ^90% and ^30% of the total continuum. In the most 
polarized sources, we also measure the broad line emission 
from the quasar with approximately the expected equiva- 
lent and velocity widths. 

- A nebular continuum which, based on the observed strength 
of the Hell A1640 recombination line, is typically ~+0% of 
the total flux in the UV. 

- The remaining continuum is assumed to arise from a young 
stellar population although we have only very weak di- 
rect observational evidence that this is the case. Constraints 
from the observed continuum colour severely restrict the 
age and reddening of such a population in the low polariza- 
tion objects. 

- Strong, spatially-extended narrow emission lines from the 
ISM of the host galaxy. From photoionization modelling 
studies it is shown that the dominant fraction of this line 
emission is excited by the ionizing radiation from the 
quasar. 

The clumpy dust-scattering model we have employed 
shows that a few percent of the hidden quasar light is rendered 
visible in the restframe UV spectrum of the radio galaxies. For 
luminous AGN, this is generally sufficient to dominate the ob- 
served continuum in this wavelength range. Only young, lu- 
minous and relatively unreddened starbursts will reveal them- 
selves clearly against the scattered 'white haze'. Given the 
large uncertainty in the reddening of starbursts and the conse- 
quent broad range in the allowed conversion of stellar UV flux 
to SFR, it is likely that we will have to await more sensitive 
sub-mm measurements in order to investigate the consistency 
between the building of the stellar content of the host spheroid 
and the chemical evolution of the ISM as inferred from the 
emission lines. 

The scattering models allow us to calculate directly the 
FIR radiation emitted from the same dust within the ioniza- 
tion cones which reflects the AGN. This is very much a lower- 
bound to the expected FIR flux since it takes account neither of 
the warm dust emission from the obscuring torus nor the cool 
dust emission from the shadowed dust outside the cones which 
could be heated by obscured star-formation as well as reflected 
quasar light. We show that the contribution from the scatter- 
ing dust is typically 10% of the observed upper limits (and one 
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850^m detection for 4C+48.48) on sub-mm emission from ra- 
dio galaxies. 

With the exception of Lya and Nv, which show substantial 
relative strength variation within our sample, the strong emis- 
sion line spectra are very similar and indicate that the ISM is 
ionized predominantly by the quasar radiation field with an ion- 
ization parameter U ~ 0.04. The presence of absorption com- 
ponents suggests that the variations in integrated Lya emission 
strength are the result of self-absorption by extended neutral 
hydrogen halos. The variation by more than a factor of 5 in 
our sample of the Nv/Civ and Nv/Hell emission line ratios 
cannot be readily explained by variations in ionization state or 
by the effects of carbon depletion onto dust. However, galac- 
tic chemical evolution models for elliptical galaxies, where the 
nitrogen enrichment arises predominantly from secondary pro- 
cesses over a timescale of ~1 Gyr, offer a good fit to the ob- 
servations of these lines and are generally consistent with the 
measurements of other lines. This scenario is similar to that 
inferred from the BLR measurements of QSO at similar red- 
shifts which imply a range of supersolar metallicities from a 
few to perhaps twenty times solar. The sequence for the radio 
galaxy ISM suggests, however, a range of metallicity extend- 
ing from about half to a few times solar. In contrast to the QSO 
BLR results, where there appear to be no metal-poor luminous 
AGN, the much more extended gas in the ISM of radio galax- 
ies shows that massive black holes can exist in relatively metal 
poor hosts. 

The Nv/Civ and Nv/Hell line ratios are positively corre- 
lated with the continuum polarization P while Lya/Civ is an- 
ticorrelated with P. We believe that this is most likely to be the 
result of an evolutionary sequence within which the spheroids 
become increasingly dust-enshrouded as the chemical enrich- 
ment of the ISM proceeds on a timescate of around 1 Gyr. A 
property of this scenario is that ULIRG like FSC 10214+4724 
are the end point rather than the beginning of such a period. 

Our relatively low resolution spectra are not very sensitive 
to the presence of weak, narrow absorption lines and so we 
cannot place very stringent constraints on the presence of lines 
arising in the photospheres of OB stars. We do, however, see 
a number of low ionization resonance and excited fine struc- 
ture lines from carbon, oxygen and silicon. These lines are typ- 
ically much weaker in absorption than seen in Lyman-break 
galaxies like MS1512-cB58 (Pettini et al. 2000 ) and, indeed, 
sometimes show redshifted emission components (see Fig. 3). 
There is a weak tendency for high polarization objects to show 
these lines in emission and low polarization objects to show 
them in absorption. This could indicate the different locations 
of the background stellar and scattered continuum sources. 

While we demonstrate that scattered quasar light generally 
dominates the UV continuum of these powerful radio galaxies, 
the presence of an evolved stellar population would be expected 
to reveal itself at wavelengths above the 4000 A break. To this 
end, we are currently observing a partially overlapping sample 
of HzRG with 2.2 < z < 2.6 in the near infrared with the 
ESO VLT in order to sample their restframe optical spectra. 
Measurements of the continuum should reveal the presence of 



any evolved stellar population while the familiar optical forbid- 
den emission lines will improve our ability to perform reliable 
abundance analyses of the ionized ISM. 
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Table 6. Emission line flux measurements. Fluxes and errors are in units of 10~ 17 ergs s _1 cm~ 2 . Excited fine structure lines are labeled with a *. Vacuum wavelengths 
are given in the second column. For blends of more than two lines, the average wavelength is shown in italics. Quoted errors for strong lines were computed as described 
in section 3.2. The typical error for weaker lines is about 10% 



□ 



Line 


\vac (A) 


4C+03.24 


0943-242 


0828+193 


4C+23.56a 


0731+438 


4C-00.54 


4C+48.48 


0211-122 


4C+40.36 


Sill 


1260.4 








3.7 




3.0 




7.5 


1.7 


Sill* 


1264.7,1265.0 




-13.4 






-4.9 


-4.6 


-12.5 






Oi, Oi*, Sra 


1304 








4.2 


3.9 


5.5 


2.1 






Sin* 


1309.3 




-10.8 


-6.8 




-8.6 






-5.7 


-9.6 


Cn , Cn* 


1334.5,1335.7 


-31 


-12.4 


-6.8 


6.7 


-9.9 


-6.1 


-6.9 


2.2 


-20.9 


Sin 


1526.7 






5.3 


3.1 






4.7 


6.8 


10.2 


Sin* 


1533.4 




-7.8 






-4.2 




2.7 


1.0 




Sin 


1808.0 




-7.9 


-18.7 




-8.4 










Sin* 


1816.9,1817.5 






-4.5 












11.6 



Table 7. Observed frame equivalent width (W\) measurements of low ionization zero-volt or excited fine-structure (labeled with a *) lines in A. W\ is positive for 
absorption lines and negative for emission lines. Vacuum wavelengths are given in the second column. For blends of more than two lines, the average wavelength is 
shown in italics. The equivalent width measurement error is typically about 10 % 
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